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ABSTRACT 

We examine the composition of matter as it flows away from gamma ray burst 
accretion disks, in order to determine what sort of nucleosynthesis may occur. 
Since there is a large flux of neutrinos leaving the surface of the disk, the electron 
fraction of the outflowing material will change due to charged current neutrino 
interactions. We calculate the electron fraction in the wind using detailed neu- 
trino fluxes from every point on the disk and study a range of trajectories and 
outflow conditions for several different accretion disk models. We find that low 
electron fractions, conducive to making r-process elements, only appear in out- 
flows from disks with high accretion rates that have a significant region both of 
trapped neutrinos and antineutrinos. Disks with lower accretion rates that have 
only a significant region of trapped neutrinos can have outflows with very high 
electron fractions, whereas the lowest accretion rate disks with little trapping 
have outflow electrons fractions of closer to one half. 

Subject headings: gamma ray:bursts-nucleosynthesis-accretion disks 

1. Introduction 

Ever since the first gamma ray bursts were detected thirty years ago, their origin has 
been a subject of great interest. Many more observations have occurred in recent years 
which point to exotic supernova as their astrophysical source; for a review see Meszaros 
(2002). Although the hydrodynamic details of the evolution of these objects is still under 
development, see e.g. Woosley (1993); MacFadyen & Woosley (1999), it is likely that the 
burst will originate from a configuration where an accretion disk surrounds a black hole. 
Several models for these accretion disks have been examined (Popham, Woosley, & Fryer 
1999; Narayan, Piran, & Kumar 2001; DiMatteo, Perna, & Narayan 2002). In addition 
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to the highly coUimated ultrarelativistic jet that produces the observed gamma rays, there 
will also be outflow from the accretion disk which has been studied by MacFadyen (2003); 
MacFadyen & Woosley (1999). This material begins at relatively high temperature (a few 
MeV) and therefore will undergo a primary nucleosynthesis where free nucleons form heavy 
nuclei. Since this is a different environment than that which occurs in the relativistic jet, 
different nucleosynthesis products will be produced. 

Any analysis of the nucleosynthesis must begin with the evolution of the electron fraction 
in the accretion disk. This was done in Pruet, Woosley, & Hoffman (2003) by considering 
only electron and positron capture and by Surman & McLaughlin (2004) by considering 
neutrino and antineutrino capture as well. The disk is hot enough that for high accretion 
rates the neutrinos and even the antineutrinos can become trapped, creating a neutrino torus 
around the black hole that is in some ways similar to the neutrinosphere at the surface of 
the protoneutron star in a normal supernova. 

A prehminary analysis of the nucleosynthesis from the outflow of accretion disks was 
done in Pruet, Thompson, & Hoffman (2004) using a spherically symmetric neutrino driven 
wind model, and by Pruet, Surman, & McLaughlin (2004) using an outflow at flxed veloc- 
ity and entropy. Maeda & Nomoto (2003) considered the nucleosynthesis that will occur 
from explosive burning in gamma ray bursts. Pruet, Thompson, & Hoffman (2004) con- 
cluded that iron peak elements will be produced in this outflow, while Pruet, Surman, & 
McLaughlin (2004) examined various rare isotopes that would point to GRB disk outflow as 
a unique nucleosynthesis event. Nucleosynthesis from accretion disks has also been discussed 
in Fujimoto et al. (2004). 

In this paper we examine the nucleosynthesis in the gamma ray burst cjccta by studying 
the effect of the neutrino ffux on the electron fraction of the material as it leaves the surface 
of the disk. We examine the conditions under which the electron fraction is quite low and the 
heaviest, r-proccss elements are likely to be produced. We also examine the the conditions 
under which the lighter iron peak nuclei are likely to be formed. 

2. Trajectories 

Determining the trajectories followed by mass elements ejected from the gamma ray 
burst accretion disk is a complicated problem. The trajectories depend on the initial ther- 
modynamic parameters of the material, and the manner in which energy is imparted to the 
gas. With different accretion disk parameters, such as the viscosity, a, spin parameter a and 
size of the black hole, the results are hkely to vary. Although disk winds have been studied 
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extensively in the context of other astrophysical objects such as AGNs and X-ray stars, a 
complete discussion of the outflow from gamma ray burst accretion disks has not yet been 
attempted. 

In order to investigate the impact of neutrino interactions on the electron fraction and 
therefore the nucleosynthesis coming from gamma ray burst disk outflow, we approximate 
disk outflow trajectories by using a two component approach. At very large distance, we use 
spherical symmetry. Initially, however, the forces pushing the material off the accretion disk 
will be closer to vertical and we assume cylindrical symmetry. At short distances, therefore, 
the material moves approximately vertically, whereas at long distances, it moves radially 
away from the black hole. 

The general hydrodynamic equations in steady state for matter conservation, momentum 
conservation and energy conservation are, e.g. Duncan, Shapiro, & Wasserman (1986): 

| = -V(PU) (1, 
Du 1 

- = --yP-y^ (2) 

i^--^-^(V.S...V.S„) (3) 

Here p is the density, P is pressure, u is the velocity, e is the energy per unit mass and 
Sph and Sj, are the photon and neutrino fluxes respectively. The gravitational potential is 
represented by 0, while the notation D/Dt represents a convective derivative. 

As discussed above, we divide the outflow into two regions, and approximate the outflow 
in these regions as having a cylindrical or spherical symmetry. For the case of spherically 
radial flow, we rewrite the mass and the momentum equation, so we have: 

M = Anr'^pur (4) 

du IdP GM 

= — ^ Y- ^ 

or p or r^ 

In these equations M is the mass loss rate, G is the gravitational constant, and M is 
the mass enclosed within radius r. We also rewrite these equations for cylindrical geometry 
for vertical flow off the disk: 



M = -27r 



pUj-rcdz + / puzrcdr, 



(6) 
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In the above equations Tc is the radial cyhndrical coordinate. 

In principle one should solve Eqs. 6-8 (or Eqs. 4-5) together with with Eq. 3. 
However, because of the complex geometry, this involves a lengthy numerical calculation 
which is not practical until more accurate disk models become available. Still, we would 
like to understand the importance of neutrino interactions in the outflow and what sort of 
electron fractions can obtain in this environment. Therefore, in place of Eqs. 3 and Eq. 7 
-8 (or Eq. 5), we use a parameterization for the velocity. 



where R = {z'^ + r'^)^'^ for the first, vertical part of the trajectory and the starting position of 
the material is Rq. Once we switch to spherical flow R = r. We study the results as functions 
of the parameters f3 = 0.2 to 2.5 , foo = 5 x lO'^kms"^ — 5 x lO'^kms^^, and Rq = 50 km to 
600 km. Although all trajectories with the same v^o asymptote to the same value at large 
distance, the ones with smaller (3 have a greater initial acceleration and arrive there more 
quickly. 

Three velocity trajectories arc shown in Fig. 1 as a function of distance from the black 
hole, all for the same p and final velocity Voo = 3 x lO^kms^^, but for two different initial 
starting points. 

We use the mass conservation equations (Eq. 8 or Eq. 5) in order to determine the 
density given the position and velocity. In the cylindrical case we assume the velocity is 
completely in the z direction, and therefore there is no expansion of the material in the 
radial direction along the disk. 

In our calculations we assume a constant entropy, since any heating of the material 
should be done at the surface of the disk. We take the entropy as an input parameter in the 
range of s = 10 to s = 40, from which we calculate the temperature at each step using the 
expression for the entropy in units of the Boltzmann constant 



where is the photon entropy, Sg+e- is the entropy of the electron-positron pairs, and Snuci 




(9) 
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is the entropy of the nucleons. The photon entropy is given by 

= 0.019^^ (11) 
Pio 

where TmcV is the temperature in units of MeV and pio is the density in units of 10^° g cm~^. 
The entropy of the electron positron pairs is 



^ 0.0022^^ 

Pio 



TmcV J \ TmcV 



(12) 



, TMeV J Jo 1 + exp(y - fjLe/TMev) 

where 7j = {x"^ + {rrie/TMev)'^)^'^, ""^e is the mass of the electron and fi^ is the chemical 
potential of the electrons, which is determined by the temperature, density and electron 
fraction of the material 



PioYe = 2.2 X 10 

where 



2m ~ ^"im 



2m 



MeV / \ MeV 



(14) 



(15) 



TueV J Jo 1 + CXp(?/ - He/TMev) ' 

We take the complete expression for the electron and positron pairs since we are neither in the 
limit where they are fully relativistic or fully non-relativistic. The approximate expression 
for the entropy of the nucleons is 

Weon = 7.4 + ln(^j. (16) 

We begin the calculations when the material is at the surface of the disk. We start 
with disk conditions from the disk models of DiMatteo, Perna, & Narayan (2002) (hereafter 
DPN) for disks with accretion rates M > 1 Mq/s and from Popham, Woosley, & Fryer (1999) 
for more slowly accreting disks. We take the disk surface to be at the density scale height 

H-\-^\-\ (17) 
paz 

We begin the outflow in the vertical direction so Ur is zero, and = \u\. We take steps in 
vertical distance, at each point determining a velocity, a density, temperature and electron 
chemical potential. We turn from the vertical solution to the cylindrical one when the 
material has reached one, two or three vertical scale heights above the disk, although this is 
the least sensitive parameter in this model. 
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The results of this formulation for a relatively high accretion rate disk model and wind 

parameters, f3 = 0.8 and = 100km, s = 20 (Model 1), = 250 km, s = 10 (Model 2) and 
Tc = 250 km, s = 20 (Model 3) are shown in Figs. 2, 3 and 4. In Fig. 2, we plot two sample 
trajectories in position space. Since we have fixed the turnover point at two vertical scale 
heights, the final trajectories appear to lie almost on top of each other. As mentioned above, 
we found little impact on the results when examining different turnover points. In Fig. 3, 
we show the densities for these trajectories. In Fig. 4 the effect of a larger entropy can also 
be seen as a greater temperature at a given distance. 



3. Neutrino Fluxes 

Next we use the trajectories developed in the previous section to calculate the evolution 
of the electron fraction due to electron, position, neutrino and antineutrino annihilation: 

e"" + p <-> i^e + n (18) 

e'^ + n ^ Ue + P (19) 

Since the electrons and positrons are in equilibrium with the baryons and photons and 
therefore their distribution can be described with the temperature, T, and electron chemical 
potential fie, the forward rates in Eqs.l8 and 19 are easy to calculate. However, the neutrinos 
are not in equilibrium and their fiux at each point on the trajectory must be calculated by 
summing over the neutrinos which originate at every part of the disk. 

In Surman & McLaughlin (2004), we calculated the evolution of the electron fraction 
of a mass element as it spiraled from the outer edge of the toward the center. This was 
done by summing all contributions form the neutrino flux at all points on the disk, taking 
into account regions where the neutrinos become optically thick. These calculated electron 
fractions are the starting Y^s for the material ejected from the disk. Furthermore, these 
neutrino fluxes are used to calculate the spectrum at every point above the disk through 
which the ejected material passes. 

In Fig. 5 we show the results of using the reactions Eqs. 18 and 19 to determine 
the electron fraction for Models 1, 2, and 3. Dashed lines show the true electron fraction 
while the dot-dashed lines show the electron fraction calculated without neutrino capture 
interactions. 

In the high entropy models (s = 20) the large electron fractions are due to the increased 
importance of electron and positron capture at higher temperature. With higher tempera- 
tures the electrons and positrons essentially reset the electron fraction to a new equilibrium 
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value and the material "forgets" its original disk value and winds up consisting of nearly 
equal numbers of neutrons and protons. This is evident in the initial sharp increase in the 
electron fraction. 

However, when the neutrino interactions have been included, the situation is quite dif- 
ferent. Even at these high entropies, the neutrinos have a marked impact on the evolution 
of the electron fraction, bringing it down to < 0.3. The reverse rates in Eq. 18 and 19 
overwhelm the forward rates. Furthermore, due to the higher energies of the electron an- 
tineutrinos, the antineutrino capture on protons is larger than neutrino capture on neutrons. 
As can be seen from the figure, the release point of the material determines the degree to 
which neutrino capture influences the final electron fractions. 

For lower entropies we find less of an increase in the electron fraction when the neutrino 
capture reactions are turned off. The lower entropy means that the system is partially 
electron degenerate, and so even while the electron and positron captures are maintaining a 
weak equilibrium by themselves, the electron fraction is quite low. Note that the low entropy 
model corresponds to essentially no heating of the outflow material since the disk has an 
entropy of order 10. 

4. Results 

In this section we explore the qualitative effect of variables such as entropy and outffow 
timescale on the electron fraction and therefore on the nucleosynthesis. We also explore 
different disk models. 

In the previous section we discussed the importance of entropy on the electron fraction. 
The outflowing material becomes less neutron rich if it is exposed to more positron capture, 
which happens when the entropy rises and the chemical potential decreases. Fig. 6 shows 
the electron fraction measured at Tg = 10 for several different entropies, for trajectories 
that start at Tq = 250 km and two different outflow parameters /3. This flgure shows that at 
the highest entropies s = 40 and fast accelerations (/? = 0.8), where the neutrinos have the 
least influence, the electron fraction can become as high as 0.45, while for slow accelerations 
(/? = 2.5) the neutrinos insure that the electron fraction is very low (~ 0.1). 

The effect of the time scale of the outflow through the parameter /3 is shown in Fig. 7. 
The material accelerates much more quickly at lower providing less opportunity for the 
neutrinos to move the system toward weak equilibrium. In the case of the high entropy, the 
effect is most pronounced. The higher betas mean more time for the neutrinos to drive the 
electron fraction down. In the case of the lower entropy, as previously noted the electrons 
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and positrons already favor a low electron fraction, and so the electron fraction never gets 
very high. 

So far we have considered disk models with accretion rates of M = IOMq/s, and spin 
parameter a — 0. This model is desirable because it produces many neutrinos which create 
a large energy deposition from neutrino-antineutrino annihilation. Such high accretion rates 
may be expected in the case of neutron star-ncutron star mergers, but lower accretion rates 
are suggested for the collapsar model (MacFadyen & Woosley 1999). However, increasing 
the spin parameter or viscosity is similar to increasing the accretion rate as far as driving up 
the neutrino flux. So the high accretion rate model we have considered may mimic a lower 
accretion rate model with a large spin parameter. 

When the accretion rate becomes lower an interesting effect occurs. At M = 10 Mq/ s, 
a = the neutrino surface is at 200 km but by M = IMq/ s, a = 0, it has shrunk to 40 
km. Similarly the antineutrino sphere shrinks from 140 km to around 32 km; see figures 3 
and 4 in Surman & McLaughhn (2004). There are relatively few antineutrinos because the 
antineutrino surface is quite small, however the neutrino capture rates are still large and the 
net effect can be to drive to very high values. This can be seen in Fig. 8 which gives 
the electron fraction for a M = 1 M©/ s model for various entropies both with neutrinos and 
without neutrinos. Note the difference in the slow acceleration (/3 = 2.5) curves in Figs. 6 
and 8. 

This model is most sensitive to the outflow parameters because the neutrino capture 
rates, while still quite strong, are smaller than in the M — 10 M©/ s model, by a factor of 
2-3 for the neutrinos and an order of magnitude for the antineutrinos. Because of the 
extreme sensitively to the neutrino parameters a neutrino diffusion calculation is needed to 
better determine the spectra and luminosity of the neutrinos and antineutrinos emitted at 
each part of the accretion disk. 

In Fig. 9, electron fractions for models with lower accretion rate disks of M = 0.1 Mq/ s, 
a = 0.95 are shown. In this model the material starts off with a low density and temperature 
relative to that of the higher accretion rate models. Although there is only a very small region 
where the neutrinos are trapped in this disk, and the antineutrinos are not trapped at all, 
the influence of neutrino capture is apparent. Even here, the electron fraction changes by as 
much as 30% in the upward direction depending on the model. Such changes will have an 
important impact on the nucleosynthesis in the iron peak region. 

Since the parameterization discussed here is independent of the heating mechanism, 
these results can be used as a rough gauge of the electron fraction and therefore the nucle- 
osynthesis for any type of wind model. 
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4.1. Preliminciry Nucleosynthesis Calculation 

The nucleosynthesis that may be produced from configurations that give a low electron 
fraction is shown in figure 10. Here we show the results of an r-process calculation for s = 10, 
(3 — 0.8 and Tq = 250 km. For this choice of parameters, the outfiow does result in r-process 
abundances. However, it can be seen from the figure that although the neutrinos initially 
helped to keep the electron fraction low, they actually cause the A=195 peak to disappear 
in an alpha effect (Fuller & Meyer 1995; McLaughhn, Fuller, & Wilson 1996). 

The neutrino driven wind models for low accretion rate disks (M = O.IMq /s) were 
discussed in Pruet, Thompson, & Hoffman (2004). These fall at around P — 2.5, v^o — 
3 X lO^kms"^ in our parameterization and we estimate that the electron fraction may be 
increased by as much as 5% - 20% by the neutrinos. These conditions may lead to a large 
overproduction of elements such as ^^Ca and ^^Sc, ^^Ti, ^^Ti, ^^Cu, ^^Zn as discussed in 
Pruet, Surman, & McLaughlin (2004). 

For winds that accelerate quickly from an accretion disk of around M — 1Mq/s and 
a = 0, the electron fraction can become as high as 0.8. In this case the nucleosynthesis would 
be dominated by nickel, in addition to making nuclei on the proton rich side of the valley of 
beta stability, such as ^®Cu, ^^Zn,^°Fe and ^^Fe. 

5. Conclusions 

We have presented a parameter study for outflow from gamma ray burst accretion 
disk and calculated the electron fractions produced in these outflows. We have considered 
the impact of charge changing neutrino interactions on the outfiow, using our previously 
calculated neutrino fiuxes from every point on the disk (Surman & McLaughlin 2004). 

Complete hydrodynamic models for the outflow for various disk parameters will become 
available in the future and the parameter study presented here can then be used to determine 
which disk models and outflows are likely to produce different types of nucleosynthesis, 
such as r-process or iron peak nuclei. For example, we find that the conditions that are 
most conducive to making the r-process elements come from those disks with high accretion 
rates or spin parameters, e.g. M = IOMq/s, a = 0, such that they produce a sizable 
region of trapped antineutrinos. Disk models with somewhat lower accretion rates and spin 
parameters, e. g. M = 1Mq/s, a — may still have a significant region of trapped neutrinos, 
but a smaller region of trapped antineutrinos. This causes the electron fractions to become 
very high, potentially as high as 0.8. Still lower accretion rate models, e.g. M — O.IMq/s, 
a = 0.95 will have have very small region of trapped neutrinos, which can raise the electron 
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fr action by 5% to 30%. 

Several parameters determine the electron fraction in the outflow, the outflow timescale, 
the entropy and the release point on the disk. The effect of the high entropy is to raise the 
electron fraction through positron capture. A more slowly accelerating wind increases the 
importance of neutrino and antineutrino capture which can drive the electron fraction up or 
down. Releasing the material closer to the center of the disk has a similar effect. 

In all cases there is interesting nucleosynthesis to explore. However, if the r-process is to 
come from gamma ray bursts, it is necessary for some of the outflow material to have a small 
electron fraction. This requires high accretion rate, high spin models, and small amounts of 
heating in the wind. 

Future observations of emission lines from gamma ray bursts, when combined with 
studies such as this one may be an avenue toward understanding not only the nucleosynthesis 
originating from GRBs but also the conditions in the outflow and in the accretion disk itself. 
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Table 1. Disk and Wind Parameters for Three Sample Trajectories 
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Fig. 1. — Shows velocity plotted against radius for three different trajectories. The sohd hne 
shows Model 1 which starts at a distance on the disk from the center of r = 100 km. The 
dashed hne shows Model 2 and 3 which start at a distance of r = 250 km from the center. 
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Fig. 2. — Shows vertical height as a function of radial coordinate, r^, for the same three 
models as in Fig. 1. The dotted line shows the scale height of the disk. 
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Fig. 3. — Shows density as a function of distance from the center, R = 
same three models as in Fig 1. 



{z^ + r^)°-^, for the 
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Fig. 4. — Shows temperature as a function of R, for Model 1 (solid line), Model 2 (short 
dashed hne) and Model 3 (long dashed hne). 
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Fig. 5. — Shows Ye as a function of distance from the black hole for the same three models as 
in Fig. 1. The dotted and dot-dashed lines show the effect when neutrino capture interactions 
are turned off. 
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Fig. 6. — For disk models with M = 10 Mq s~^, a = 0. Shows Y^. as a function of s, using a 
release point of rg = 250 km and a final velocity of Voo = 3 x lO'^kms"^. The neutrinos have 
the maximum effect in a slowly accelerating outflow (/3 = 2.5) and significantly decrease the 
electron fraction. 
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Fig. 8. — For a disk model with an accretion rate of M = 1.0 M0 s^-*^, a = 0. Shows as a 
function of s, using a release point of = 170km and a final velocity v^o = 3 x lO^kms"^. 
The neutrinos raise the electron fraction considerably in a slowly accelerating outflow {/3 — 
2.5). 
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Fig. 9. — For a disk model with an accretion rate of M = 0.1 M0 s~^, spin parameter 
a = 0.95. Shows Yg as a function of s, using a release point of = 100 km and a final 
velocity Voo = 3 x lO^km s~^ . The neutrinos raise the electron fraction. For slowly accelerating 
outflows (/3 = 2.5), the effect can be as large as 30%. 
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Fig. 10. — For a disk model with an accretion rate of M = 10 Mq s^^, spin parameter a — 
and for a trajectory of s = 10, /3 = 0.8, v^o = 1 x lO^'kms^^ and tq = 250 km. The sohd 
hne gives the result with all neutrino interactions in the wind turned on, the dashed line has 
neutrino interactions in the disk, but not in the wind. 



